Scripts {pliman} 67ª RBRAS e 20º SEAGRO

1 Pacotes

library(pliman)
library(tidyverse)
setwd("D:/Desktop/UFSC/cursos/pliman67rbras")
set_wd_here()
# set_pliman_viewer("mapview")

2 Importação e manipulação

2.1 Importação

leaves <- image_import("leaves.jpg", plot = TRUE)

2.2 Resolução

dpi(leaves, viewer = "mapview")

2.3 Cortar, girar e redimensionar

crop <- image_autocrop(leaves, plot = TRUE)

# crop <- image_crop(leaves, viewer = "mapview")

crop <- image_crop(leaves,
                   plot = TRUE)
## Error in round(minw, 0):round(maxw, 0): resultado seria um vetor muito longo

3 Segmentação

# Índices para segmentação
pliman_indexes()
##  [1] "R"     "G"     "B"     "NR"    "NG"    "NB"    "GB"    "RB"    "GR"   
## [10] "BI"    "BIM"   "SCI"   "GLI"   "HI"    "NGRDI" "NDGBI" "NDRBI" "I"    
## [19] "S"     "VARI"  "HUE"   "HUE2"  "BGI"   "L"     "GRAY"  "GLAI"  "CI"   
## [28] "SHP"   "RI"    "G-B"   "G-R"   "R-G"   "R-B"   "B-R"   "B-G"   "DGCI" 
## [37] "GRAY2" "L*"    "a"     "b*"    "L*-a"  "L*-b"  "b*-a"
index <- image_index(leaves, index = c("R", "G", "B", "B-R"))

plot(index, type = "density")

# Imagem binária
binary <- image_binary(leaves, index = "B-R")

binary <- 
  image_binary(leaves,
               index = "B-R",
               invert = TRUE)

binary <- 
  image_binary(leaves, 
               index = "B-R",
               fill_hull = TRUE)

binary <- 
  image_binary(leaves, 
               index = "B-R",
               filter = 5,
               fill_hull = TRUE)

# Segmentação
seg <- 
  image_segment(leaves,
                # threshold = 0.1,
                # threshold = "aa",
                index = "B-R")

# Threshold adaptativo
adaptive <- image_import("adaptive.jpg", plot = TRUE)

# Ìndice de imagem
ind <- image_index(adaptive, index = "R")

plot(ind, type = "density")

# Segmentação por Otsu
image_segment(adaptive, index = "R")

# Segmentação utilizando threshold adaptativo
image_segment(adaptive,
              index = "R",
              # filter = 5,
              # fill_hull = TRUE,
              threshold = "adaptive")

4 Análise de objetos

4.1 Imagem única

grains <- image_import("grains.jpg", plot = TRUE)

# Analisar os objetos
res <- analyze_objects(grains, index = "GRAY")

res$statistics
# Algumas funcionalidades
res <- 
  analyze_objects(grains, 
                  index = "GRAY",
                  lower_noise = 0.4,      # remove ruídos maiores
                  show_contour = FALSE,
                  marker = "point")

plot(res)

res <- 
  analyze_objects(grains, 
                  index = "GRAY",
                  topn_upper = 10,        # Somente os 10 maiores grãos
                  show_contour = FALSE,
                  marker = "point")

4.2 Processamento em lote com correção (referência)

4.2.1 Pelo tamanho da referência

# Analisar os objetos
res <- 
  analyze_objects(pattern = "P",
                  dir_original = "flax_grain",    # subpasta com as imagens originais
                  dir_processed = "proc",         # subpasta com as imagens processadas
                  reference = TRUE,               # indica que há uma referência
                  reference_larger = TRUE,        # indica que a referência é o maior objeto 
                  reference_area = 6,             # a referência tem 6 cm2
                  index = "GRAY",
                  show_contour = FALSE,
                  marker = "point",
                  plot = FALSE,
                  save_image = TRUE)
## Processing image P01 |=======                                    | 17% 00:00:00 
## Processing image P02 |==============                             | 33% 00:00:03 
## Processing image P03 |======================                     | 50% 00:00:07 
## Processing image P04 |=============================              | 67% 00:00:11 
## Processing image P05 |====================================       | 83% 00:00:13 
## Processing image P06 |===========================================| 100% 00:00:16 
## --------------------------------------------
##  Image Objects
##    P01     277
##    P02      91
##    P03     179
##    P04     112
##    P05      22
##    P06      78
## --------------------------------------------
hist(res$results$length)

4.2.2 Pela cor da referência

flax <- image_import("flax_af/A2_32_3.jpg", plot = TRUE)

# Índice para segmentar a referência e folhas do fundo
image_segment_iter(flax, 
                   index = c("R/(G/B)", "B-R"),
                   ncol = 3)
##      image pixels   percent
## 1 original 806490 100.00000
## 2     seg1 129657  16.07670
## 3     seg2  38517  29.70684

# Note que agora o processamento é realizado de forma paralela
res <-
  analyze_objects(pattern = "A",
                  dir_original = "flax_af",
                  reference = TRUE,
                  reference_area = 20,
                  back_fore_index = "R/(G/B)",
                  fore_ref_index = "B-R",
                  watershed = FALSE,
                  filter = 2,
                  plot = FALSE,
                  parallel = TRUE)
## --------------------------------------------
##      Image Objects
##    A1_28_1      18
##    A1_28_2      28
##    A1_28_3      28
##   A10_90_1     149
##   A10_90_2     116
##   A10_90_3      77
##   A11_98_1     132
##   A11_98_2     151
##   A11_98_3     119
##  A12_105_1     110
##  A12_105_2     127
##  A12_105_3     154
##    A2_32_1      23
##    A2_32_2      25
##    A2_32_3      31
##    A3_42_1      39
##    A3_42_2      47
##    A3_42_3      52
##    A4_46_1      58
##    A4_46_2      38
##    A4_46_3      56
##    A5_55_1      53
##    A5_55_2      74
##    A5_55_3      61
##    A6_63_1      68
##    A6_63_2     105
##    A6_63_3      83
##    A7_70_1      98
##    A7_70_2     111
##    A7_70_3     112
##    A8_76_1     113
##    A8_76_2     118
##    A8_76_3     114
##    A9_83_1     120
##    A9_83_2      89
##    A9_83_3     101
## --------------------------------------------
merged <- get_measures(res)


# Organizar os dados para ajuste do modelo Logístico
df_plot <- 
  merged$summary |> 
  separate_col(img, 
               into = c("avaliacao", "das", "bloco")) |> 
  mutate(das = as.numeric(das))

formula <- y ~ b1/(1 + exp(b2 - b3 * x))

ggplot(df_plot, aes(das, area_sum)) + 
  geom_smooth(method = "nls",
              method.args = list(formula = formula,
                                 start = c(b1 = 248,
                                           b2 = 6,
                                           b3 = 0.07)),
              se = FALSE,
              color = "red") +
  stat_summary(fun.data = mean_se,
               geom = "errorbar",
               width = 0.5) +
  stat_summary(fun = mean,
               geom = "point",
               col = "blue",
               size = 3) +
  scale_x_continuous(breaks = unique(df_plot$das)) +
  scale_y_continuous(breaks = seq(0, 150, by = 25)) +
  labs(x = "Dias após a semeadura",
       y = expression(Área~foliar~média~(cm^2~planta^{-1}))) +
  theme_bw(base_size = 16) +
  theme(panel.grid = element_blank()) 

5 Ortomosaicos

Imagens disponibilizadas por Filipe Matias, no github do FIELDimageR

5.1 Cobertura vegetal (canopy)

# imagem
img <- 
  image_import("potatoes.tif", path = "orthomosaic") |> 
  image_horizontal()

plot(img)

# preparar a imagem
# Somente iterativo
# prep <- image_prepare_mv(img)
prep <- 
  img |> 
  image_rotate(-2.076, plot = FALSE) |> 
  image_crop(width = 194:1075,
             height = 143:379,
             plot = TRUE)

# índice para segmentação
image_index(prep, index =  "HUE")

res <- 
  analyze_objects_shp(prep,
                      nrow = 16,
                      ncol = 9,
                      index = "HUE",
                      object_index = c("NGRDI", "DGCI"))

# cobertura de solo
plot_index_shp(res)

# resumir as informações
results <- get_measures(res)

5.2 Índices de vegetação

# Nível de parcela
plot_index_shp(res, attribute = "NGRDI")

# Nível de pixel
plot_index(object = res, index = "DGCI")

# plot_index(object = res, index = "DGCI", viewer = "map")

# resumir as informações
results <- get_measures(res)

5.3 Stand de plantas

# imagem
stand <- image_import("stand.jpg", path = "orthomosaic", plot = TRUE)

# índice para segmentação
image_segment(stand, "HUE")

image_segment(stand, "NGRDI")

image_segment(stand, "NGRDI", invert = TRUE)

res <- 
  analyze_objects_shp(stand,
                      index = "NGRDI",
                      invert = TRUE,
                      nrow = 7,
                      ncol = 1)

par(mfrow = c(1, 2))
# cobertura de solo
plot_index_shp(res)

# mapeamento das plantas (distâncias e CV(%))
object_mark(res)
mapped <-  object_map(res)

# coeficiente de variação dentro da linha
barplot(mapped$cvs)

par(mfrow = c(1, 1))

6 Fitopatometria

6.1 Iterativa

folha <- image_pliman("sev_leaf.jpg", plot = TRUE)

# measure_disease_iter(folha, viewer = "mapview")

6.2 Usando paletas

folha <- image_import("soybean_rust/soy_1.jpg", plot = TRUE)

pals <- image_import(pattern = "soja_",
                     path = "soybean_rust",
                     plot = TRUE,
                     ncol = 3)

sev <- 
  measure_disease(folha,
                  img_healthy = pals$soja_h.png,     # paleta para folha sadia
                  img_symptoms = pals$soja_s.png,    # paleta para folha doente
                  img_background = pals$soja_b.png)  # paleta para fundo

sev$severity

6.3 Usando índices de imagem

# Índice para segmentar a folha do fundo, depois a doença da folha
image_segment_iter(folha, 
                   index = c("B", "NGRDI"),
                   ncol = 3)
##      image  pixels   percent
## 1 original 4813707 100.00000
## 2     seg1 2781900  57.79122
## 3     seg2  384237  13.81203

# IMPORTANTE: theshold vai alterar a severidade
sev2 <- 
  measure_disease(folha,
                  index_lb = "B",            # índice para segmentar a folha do fundo
                  index_dh = "NGRDI",        # índice para segmentar folha doente e sadia
                  threshold = c("Otsu", 0))  # threshold para os dois índices, respectivamente

sev2$severity
# em lote (utilizando as paletas)
sev3 <- 
  measure_disease(pattern = "soy_",
                  dir_original = "soybean_rust",    # pasta com as imagens originais
                  dir_processed = "proc_rust",      # pasta com as imagens processadas
                  save_image = TRUE,
                  plot = FALSE,
                  img_healthy = pals$soja_h.png,
                  img_symptoms = pals$soja_s.png,
                  img_background = pals$soja_b.png,
                  parallel = TRUE)
# standard area diagrams
sad(sev3, n = 6)
sev4 <- 
  measure_disease(pattern = "soy_",
                  dir_original = "soybean_rust",
                  plot = FALSE,
                  index_lb = "B",
                  index_dh = "NGRDI",
                  threshold = c("Otsu", 0),
                  parallel = TRUE)


# paletas vs índices de imagem

df <- 
  data.frame(paletas = sev3$severity$symptomatic,
             índices = sev4$severity$symptomatic)
df
ggplot(df, aes(paletas, índices)) +
  geom_point() +
  geom_abline(intercept = 0, slope = 1, color = "red") +
  theme_bw(base_size = 16) +
  xlim(0, 100) +
  ylim(0, 100) +
  coord_equal()

6.4 Em grids

folha_grid <- image_import("disease_shp.jpg", plot = TRUE)

sev5 <- 
  measure_disease_shp(folha_grid,
                      index_lb = "HUE2",
                      index_dh = "NGRDI",
                      threshold = c("Otsu", 0),
                      filter = 10,
                      ncol = 3,
                      nrow = 3,
                      plot = FALSE)
plot(folha_grid)
plot(sev5$shapefiles)
plot_measures(sev5, "symptomatic")


  1. Universidade Federal de Santa Catarina, https://olivoto.netlify.app/↩︎

LS0tDQp0aXRsZTogIlNjcmlwdHMge3BsaW1hbn0gNjfCqiBSQlJBUyBlIDIwwrogU0VBR1JPIg0KYXV0aG9yOg0KICAtIFRpYWdvIE9saXZvdG9eW1VuaXZlcnNpZGFkZSBGZWRlcmFsIGRlIFNhbnRhIENhdGFyaW5hLCBodHRwczovL29saXZvdG8ubmV0bGlmeS5hcHAvXQ0Kb3V0cHV0Og0KICBybWRmb3JtYXRzOjpyZWFkdGhlZG93bjoNCiAgICBoaWdobGlnaHQ6IHRhbmdvDQogICAgZGZfcHJpbnQ6IHBhZ2VkDQogICAgbGlnaHRib3g6IHRydWUNCiAgICBnYWxsZXJ5OiB0cnVlDQogICAgbnVtYmVyX3NlY3Rpb25zOiB0cnVlDQogICAgY29kZV9kb3dubG9hZDogdHJ1ZQ0KICAgIHNlbGZfY29udGFpbmVkOiB0cnVlDQoNCi0tLQ0KDQpgYGB7ciBzZXR1cCwgaW5jbHVkZT1GQUxTRX0NCmtuaXRyOjpvcHRzX2NodW5rJHNldChlY2hvID0gVFJVRSwNCiAgICAgICAgICAgICAgICAgICAgICBlcnJvciA9IFRSVUUsDQogICAgICAgICAgICAgICAgICAgICAgd2FybmluZyA9IEZBTFNFLA0KICAgICAgICAgICAgICAgICAgICAgIG1lc3NhZ2UgPSBGQUxTRSwNCiAgICAgICAgICAgICAgICAgICAgICBjYWNoZSA9IEZBTFNFLA0KICAgICAgICAgICAgICAgICAgICAgIGZpZy5hbGlnbiA9ICJjZW50ZXIiKQ0KDQpgYGANCg0KWyFbXShsb2dvLnBuZyldKGh0dHBzOi8vNjdyYnJhczIwc2VhZ3JvLmNvbS5ici8pDQoNCiMgUGFjb3Rlcw0KDQpgYGB7cn0NCmxpYnJhcnkocGxpbWFuKQ0KbGlicmFyeSh0aWR5dmVyc2UpDQpzZXR3ZCgiRDovRGVza3RvcC9VRlNDL2N1cnNvcy9wbGltYW42N3JicmFzIikNCnNldF93ZF9oZXJlKCkNCiMgc2V0X3BsaW1hbl92aWV3ZXIoIm1hcHZpZXciKQ0KDQpgYGANCg0KIyBJbXBvcnRhw6fDo28gZSBtYW5pcHVsYcOnw6NvDQojIyBJbXBvcnRhw6fDo28NCmBgYHtyfQ0KbGVhdmVzIDwtIGltYWdlX2ltcG9ydCgibGVhdmVzLmpwZyIsIHBsb3QgPSBUUlVFKQ0KDQpgYGANCg0KIyMgUmVzb2x1w6fDo28NCmBgYHtyIGV2YWw9RkFMU0V9DQpkcGkobGVhdmVzLCB2aWV3ZXIgPSAibWFwdmlldyIpDQpgYGANCg0KDQojIyBDb3J0YXIsIGdpcmFyIGUgcmVkaW1lbnNpb25hcg0KYGBge3J9DQpjcm9wIDwtIGltYWdlX2F1dG9jcm9wKGxlYXZlcywgcGxvdCA9IFRSVUUpDQojIGNyb3AgPC0gaW1hZ2VfY3JvcChsZWF2ZXMsIHZpZXdlciA9ICJtYXB2aWV3IikNCg0KY3JvcCA8LSBpbWFnZV9jcm9wKGxlYXZlcywNCiAgICAgICAgICAgICAgICAgICBwbG90ID0gVFJVRSkNCmBgYA0KDQoNCiMgU2VnbWVudGHDp8Ojbw0KYGBge3J9DQojIMONbmRpY2VzIHBhcmEgc2VnbWVudGHDp8Ojbw0KcGxpbWFuX2luZGV4ZXMoKQ0KaW5kZXggPC0gaW1hZ2VfaW5kZXgobGVhdmVzLCBpbmRleCA9IGMoIlIiLCAiRyIsICJCIiwgIkItUiIpKQ0KcGxvdChpbmRleCwgdHlwZSA9ICJkZW5zaXR5IikNCg0KIyBJbWFnZW0gYmluw6FyaWENCmJpbmFyeSA8LSBpbWFnZV9iaW5hcnkobGVhdmVzLCBpbmRleCA9ICJCLVIiKQ0KYmluYXJ5IDwtIA0KICBpbWFnZV9iaW5hcnkobGVhdmVzLA0KICAgICAgICAgICAgICAgaW5kZXggPSAiQi1SIiwNCiAgICAgICAgICAgICAgIGludmVydCA9IFRSVUUpDQpiaW5hcnkgPC0gDQogIGltYWdlX2JpbmFyeShsZWF2ZXMsIA0KICAgICAgICAgICAgICAgaW5kZXggPSAiQi1SIiwNCiAgICAgICAgICAgICAgIGZpbGxfaHVsbCA9IFRSVUUpDQpiaW5hcnkgPC0gDQogIGltYWdlX2JpbmFyeShsZWF2ZXMsIA0KICAgICAgICAgICAgICAgaW5kZXggPSAiQi1SIiwNCiAgICAgICAgICAgICAgIGZpbHRlciA9IDUsDQogICAgICAgICAgICAgICBmaWxsX2h1bGwgPSBUUlVFKQ0KDQojIFNlZ21lbnRhw6fDo28NCnNlZyA8LSANCiAgaW1hZ2Vfc2VnbWVudChsZWF2ZXMsDQogICAgICAgICAgICAgICAgIyB0aHJlc2hvbGQgPSAwLjEsDQogICAgICAgICAgICAgICAgIyB0aHJlc2hvbGQgPSAiYWEiLA0KICAgICAgICAgICAgICAgIGluZGV4ID0gIkItUiIpDQoNCiMgVGhyZXNob2xkIGFkYXB0YXRpdm8NCmFkYXB0aXZlIDwtIGltYWdlX2ltcG9ydCgiYWRhcHRpdmUuanBnIiwgcGxvdCA9IFRSVUUpDQoNCiMgw4xuZGljZSBkZSBpbWFnZW0NCmluZCA8LSBpbWFnZV9pbmRleChhZGFwdGl2ZSwgaW5kZXggPSAiUiIpDQpwbG90KGluZCwgdHlwZSA9ICJkZW5zaXR5IikNCg0KIyBTZWdtZW50YcOnw6NvIHBvciBPdHN1DQppbWFnZV9zZWdtZW50KGFkYXB0aXZlLCBpbmRleCA9ICJSIikNCg0KIyBTZWdtZW50YcOnw6NvIHV0aWxpemFuZG8gdGhyZXNob2xkIGFkYXB0YXRpdm8NCmltYWdlX3NlZ21lbnQoYWRhcHRpdmUsDQogICAgICAgICAgICAgIGluZGV4ID0gIlIiLA0KICAgICAgICAgICAgICAjIGZpbHRlciA9IDUsDQogICAgICAgICAgICAgICMgZmlsbF9odWxsID0gVFJVRSwNCiAgICAgICAgICAgICAgdGhyZXNob2xkID0gImFkYXB0aXZlIikNCmBgYA0KDQoNCg0KIyBBbsOhbGlzZSBkZSBvYmpldG9zDQojIyBJbWFnZW0gw7puaWNhDQpgYGB7cn0NCmdyYWlucyA8LSBpbWFnZV9pbXBvcnQoImdyYWlucy5qcGciLCBwbG90ID0gVFJVRSkNCg0KIyBBbmFsaXNhciBvcyBvYmpldG9zDQpyZXMgPC0gYW5hbHl6ZV9vYmplY3RzKGdyYWlucywgaW5kZXggPSAiR1JBWSIpDQpyZXMkc3RhdGlzdGljcw0KDQojIEFsZ3VtYXMgZnVuY2lvbmFsaWRhZGVzDQpyZXMgPC0gDQogIGFuYWx5emVfb2JqZWN0cyhncmFpbnMsIA0KICAgICAgICAgICAgICAgICAgaW5kZXggPSAiR1JBWSIsDQogICAgICAgICAgICAgICAgICBsb3dlcl9ub2lzZSA9IDAuNCwgICAgICAjIHJlbW92ZSBydcOtZG9zIG1haW9yZXMNCiAgICAgICAgICAgICAgICAgIHNob3dfY29udG91ciA9IEZBTFNFLA0KICAgICAgICAgICAgICAgICAgbWFya2VyID0gInBvaW50IikNCnBsb3QocmVzKQ0KDQpyZXMgPC0gDQogIGFuYWx5emVfb2JqZWN0cyhncmFpbnMsIA0KICAgICAgICAgICAgICAgICAgaW5kZXggPSAiR1JBWSIsDQogICAgICAgICAgICAgICAgICB0b3BuX3VwcGVyID0gMTAsICAgICAgICAjIFNvbWVudGUgb3MgMTAgbWFpb3JlcyBncsOjb3MNCiAgICAgICAgICAgICAgICAgIHNob3dfY29udG91ciA9IEZBTFNFLA0KICAgICAgICAgICAgICAgICAgbWFya2VyID0gInBvaW50IikNCg0KYGBgDQoNCg0KIyMgUHJvY2Vzc2FtZW50byBlbSBsb3RlIGNvbSBjb3JyZcOnw6NvIChyZWZlcsOqbmNpYSkNCiMjIyBQZWxvIHRhbWFuaG8gZGEgcmVmZXLDqm5jaWENCmBgYHtyfQ0KIyBBbmFsaXNhciBvcyBvYmpldG9zDQpyZXMgPC0gDQogIGFuYWx5emVfb2JqZWN0cyhwYXR0ZXJuID0gIlAiLA0KICAgICAgICAgICAgICAgICAgZGlyX29yaWdpbmFsID0gImZsYXhfZ3JhaW4iLCAgICAjIHN1YnBhc3RhIGNvbSBhcyBpbWFnZW5zIG9yaWdpbmFpcw0KICAgICAgICAgICAgICAgICAgZGlyX3Byb2Nlc3NlZCA9ICJwcm9jIiwgICAgICAgICAjIHN1YnBhc3RhIGNvbSBhcyBpbWFnZW5zIHByb2Nlc3NhZGFzDQogICAgICAgICAgICAgICAgICByZWZlcmVuY2UgPSBUUlVFLCAgICAgICAgICAgICAgICMgaW5kaWNhIHF1ZSBow6EgdW1hIHJlZmVyw6puY2lhDQogICAgICAgICAgICAgICAgICByZWZlcmVuY2VfbGFyZ2VyID0gVFJVRSwgICAgICAgICMgaW5kaWNhIHF1ZSBhIHJlZmVyw6puY2lhIMOpIG8gbWFpb3Igb2JqZXRvIA0KICAgICAgICAgICAgICAgICAgcmVmZXJlbmNlX2FyZWEgPSA2LCAgICAgICAgICAgICAjIGEgcmVmZXLDqm5jaWEgdGVtIDYgY20yDQogICAgICAgICAgICAgICAgICBpbmRleCA9ICJHUkFZIiwNCiAgICAgICAgICAgICAgICAgIHNob3dfY29udG91ciA9IEZBTFNFLA0KICAgICAgICAgICAgICAgICAgbWFya2VyID0gInBvaW50IiwNCiAgICAgICAgICAgICAgICAgIHBsb3QgPSBGQUxTRSwNCiAgICAgICAgICAgICAgICAgIHNhdmVfaW1hZ2UgPSBUUlVFKQ0KDQpoaXN0KHJlcyRyZXN1bHRzJGxlbmd0aCkNCg0KYGBgDQoNCg0KIyMjIFBlbGEgY29yIGRhIHJlZmVyw6puY2lhDQpgYGB7cn0NCiN8IG91dC13aWR0aDogIjEwMCUiDQpmbGF4IDwtIGltYWdlX2ltcG9ydCgiZmxheF9hZi9BMl8zMl8zLmpwZyIsIHBsb3QgPSBUUlVFKQ0KDQojIMONbmRpY2UgcGFyYSBzZWdtZW50YXIgYSByZWZlcsOqbmNpYSBlIGZvbGhhcyBkbyBmdW5kbw0KaW1hZ2Vfc2VnbWVudF9pdGVyKGZsYXgsIA0KICAgICAgICAgICAgICAgICAgIGluZGV4ID0gYygiUi8oRy9CKSIsICJCLVIiKSwNCiAgICAgICAgICAgICAgICAgICBuY29sID0gMykNCg0KDQojIE5vdGUgcXVlIGFnb3JhIG8gcHJvY2Vzc2FtZW50byDDqSByZWFsaXphZG8gZGUgZm9ybWEgcGFyYWxlbGENCnJlcyA8LQ0KICBhbmFseXplX29iamVjdHMocGF0dGVybiA9ICJBIiwNCiAgICAgICAgICAgICAgICAgIGRpcl9vcmlnaW5hbCA9ICJmbGF4X2FmIiwNCiAgICAgICAgICAgICAgICAgIHJlZmVyZW5jZSA9IFRSVUUsDQogICAgICAgICAgICAgICAgICByZWZlcmVuY2VfYXJlYSA9IDIwLA0KICAgICAgICAgICAgICAgICAgYmFja19mb3JlX2luZGV4ID0gIlIvKEcvQikiLA0KICAgICAgICAgICAgICAgICAgZm9yZV9yZWZfaW5kZXggPSAiQi1SIiwNCiAgICAgICAgICAgICAgICAgIHdhdGVyc2hlZCA9IEZBTFNFLA0KICAgICAgICAgICAgICAgICAgZmlsdGVyID0gMiwNCiAgICAgICAgICAgICAgICAgIHBsb3QgPSBGQUxTRSwNCiAgICAgICAgICAgICAgICAgIHBhcmFsbGVsID0gVFJVRSkNCg0KbWVyZ2VkIDwtIGdldF9tZWFzdXJlcyhyZXMpDQoNCg0KIyBPcmdhbml6YXIgb3MgZGFkb3MgcGFyYSBhanVzdGUgZG8gbW9kZWxvIExvZ8Otc3RpY28NCmRmX3Bsb3QgPC0gDQogIG1lcmdlZCRzdW1tYXJ5IHw+IA0KICBzZXBhcmF0ZV9jb2woaW1nLCANCiAgICAgICAgICAgICAgIGludG8gPSBjKCJhdmFsaWFjYW8iLCAiZGFzIiwgImJsb2NvIikpIHw+IA0KICBtdXRhdGUoZGFzID0gYXMubnVtZXJpYyhkYXMpKQ0KDQpmb3JtdWxhIDwtIHkgfiBiMS8oMSArIGV4cChiMiAtIGIzICogeCkpDQoNCmdncGxvdChkZl9wbG90LCBhZXMoZGFzLCBhcmVhX3N1bSkpICsgDQogIGdlb21fc21vb3RoKG1ldGhvZCA9ICJubHMiLA0KICAgICAgICAgICAgICBtZXRob2QuYXJncyA9IGxpc3QoZm9ybXVsYSA9IGZvcm11bGEsDQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBzdGFydCA9IGMoYjEgPSAyNDgsDQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgYjIgPSA2LA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGIzID0gMC4wNykpLA0KICAgICAgICAgICAgICBzZSA9IEZBTFNFLA0KICAgICAgICAgICAgICBjb2xvciA9ICJyZWQiKSArDQogIHN0YXRfc3VtbWFyeShmdW4uZGF0YSA9IG1lYW5fc2UsDQogICAgICAgICAgICAgICBnZW9tID0gImVycm9yYmFyIiwNCiAgICAgICAgICAgICAgIHdpZHRoID0gMC41KSArDQogIHN0YXRfc3VtbWFyeShmdW4gPSBtZWFuLA0KICAgICAgICAgICAgICAgZ2VvbSA9ICJwb2ludCIsDQogICAgICAgICAgICAgICBjb2wgPSAiYmx1ZSIsDQogICAgICAgICAgICAgICBzaXplID0gMykgKw0KICBzY2FsZV94X2NvbnRpbnVvdXMoYnJlYWtzID0gdW5pcXVlKGRmX3Bsb3QkZGFzKSkgKw0KICBzY2FsZV95X2NvbnRpbnVvdXMoYnJlYWtzID0gc2VxKDAsIDE1MCwgYnkgPSAyNSkpICsNCiAgbGFicyh4ID0gIkRpYXMgYXDDs3MgYSBzZW1lYWR1cmEiLA0KICAgICAgIHkgPSBleHByZXNzaW9uKMOBcmVhfmZvbGlhcn5tw6lkaWF+KGNtXjJ+cGxhbnRhXnstMX0pKSkgKw0KICB0aGVtZV9idyhiYXNlX3NpemUgPSAxNikgKw0KICB0aGVtZShwYW5lbC5ncmlkID0gZWxlbWVudF9ibGFuaygpKSANCg0KYGBgDQoNCg0KDQojIE9ydG9tb3NhaWNvcw0KSW1hZ2VucyBkaXNwb25pYmlsaXphZGFzIHBvciBGaWxpcGUgTWF0aWFzLCBubyBnaXRodWIgZG8gW0ZJRUxEaW1hZ2VSXShSaHR0cHM6Ly9naXRodWIuY29tL09wZW5Ecm9uZU1hcC9GSUVMRGltYWdlUikNCg0KIyMgQ29iZXJ0dXJhIHZlZ2V0YWwgKGNhbm9weSkNCg0KYGBge3J9DQojIGltYWdlbQ0KaW1nIDwtIA0KICBpbWFnZV9pbXBvcnQoInBvdGF0b2VzLnRpZiIsIHBhdGggPSAib3J0aG9tb3NhaWMiKSB8PiANCiAgaW1hZ2VfaG9yaXpvbnRhbCgpDQoNCnBsb3QoaW1nKQ0KDQojIHByZXBhcmFyIGEgaW1hZ2VtDQojIFNvbWVudGUgaXRlcmF0aXZvDQojIHByZXAgPC0gaW1hZ2VfcHJlcGFyZV9tdihpbWcpDQpwcmVwIDwtIA0KICBpbWcgfD4gDQogIGltYWdlX3JvdGF0ZSgtMi4wNzYsIHBsb3QgPSBGQUxTRSkgfD4gDQogIGltYWdlX2Nyb3Aod2lkdGggPSAxOTQ6MTA3NSwNCiAgICAgICAgICAgICBoZWlnaHQgPSAxNDM6Mzc5LA0KICAgICAgICAgICAgIHBsb3QgPSBUUlVFKQ0KDQoNCiMgw61uZGljZSBwYXJhIHNlZ21lbnRhw6fDo28NCmltYWdlX2luZGV4KHByZXAsIGluZGV4ID0gICJIVUUiKQ0KDQoNCnJlcyA8LSANCiAgYW5hbHl6ZV9vYmplY3RzX3NocChwcmVwLA0KICAgICAgICAgICAgICAgICAgICAgIG5yb3cgPSAxNiwNCiAgICAgICAgICAgICAgICAgICAgICBuY29sID0gOSwNCiAgICAgICAgICAgICAgICAgICAgICBpbmRleCA9ICJIVUUiLA0KICAgICAgICAgICAgICAgICAgICAgIG9iamVjdF9pbmRleCA9IGMoIk5HUkRJIiwgIkRHQ0kiKSkNCg0KIyBjb2JlcnR1cmEgZGUgc29sbw0KcGxvdF9pbmRleF9zaHAocmVzKQ0KDQojIHJlc3VtaXIgYXMgaW5mb3JtYcOnw7Vlcw0KcmVzdWx0cyA8LSBnZXRfbWVhc3VyZXMocmVzKQ0KYGBgDQoNCg0KDQojIyDDjW5kaWNlcyBkZSB2ZWdldGHDp8Ojbw0KDQpgYGB7cn0NCiMgTsOtdmVsIGRlIHBhcmNlbGENCnBsb3RfaW5kZXhfc2hwKHJlcywgYXR0cmlidXRlID0gIk5HUkRJIikNCg0KIyBOw612ZWwgZGUgcGl4ZWwNCnBsb3RfaW5kZXgob2JqZWN0ID0gcmVzLCBpbmRleCA9ICJER0NJIikNCiMgcGxvdF9pbmRleChvYmplY3QgPSByZXMsIGluZGV4ID0gIkRHQ0kiLCB2aWV3ZXIgPSAibWFwIikNCg0KIyByZXN1bWlyIGFzIGluZm9ybWHDp8O1ZXMNCnJlc3VsdHMgPC0gZ2V0X21lYXN1cmVzKHJlcykNCmBgYA0KDQoNCiMjIFN0YW5kIGRlIHBsYW50YXMNCg0KYGBge3J9DQojIGltYWdlbQ0Kc3RhbmQgPC0gaW1hZ2VfaW1wb3J0KCJzdGFuZC5qcGciLCBwYXRoID0gIm9ydGhvbW9zYWljIiwgcGxvdCA9IFRSVUUpDQoNCg0KIyDDrW5kaWNlIHBhcmEgc2VnbWVudGHDp8Ojbw0KaW1hZ2Vfc2VnbWVudChzdGFuZCwgIkhVRSIpDQppbWFnZV9zZWdtZW50KHN0YW5kLCAiTkdSREkiKQ0KaW1hZ2Vfc2VnbWVudChzdGFuZCwgIk5HUkRJIiwgaW52ZXJ0ID0gVFJVRSkNCg0KcmVzIDwtIA0KICBhbmFseXplX29iamVjdHNfc2hwKHN0YW5kLA0KICAgICAgICAgICAgICAgICAgICAgIGluZGV4ID0gIk5HUkRJIiwNCiAgICAgICAgICAgICAgICAgICAgICBpbnZlcnQgPSBUUlVFLA0KICAgICAgICAgICAgICAgICAgICAgIG5yb3cgPSA3LA0KICAgICAgICAgICAgICAgICAgICAgIG5jb2wgPSAxKQ0KDQpwYXIobWZyb3cgPSBjKDEsIDIpKQ0KIyBjb2JlcnR1cmEgZGUgc29sbw0KcGxvdF9pbmRleF9zaHAocmVzKQ0KDQojIG1hcGVhbWVudG8gZGFzIHBsYW50YXMgKGRpc3TDom5jaWFzIGUgQ1YoJSkpDQpvYmplY3RfbWFyayhyZXMpDQptYXBwZWQgPC0gIG9iamVjdF9tYXAocmVzKQ0KDQojIGNvZWZpY2llbnRlIGRlIHZhcmlhw6fDo28gZGVudHJvIGRhIGxpbmhhDQpiYXJwbG90KG1hcHBlZCRjdnMpDQoNCnBhcihtZnJvdyA9IGMoMSwgMSkpDQpgYGANCg0KDQoNCiMgRml0b3BhdG9tZXRyaWENCiMjIEl0ZXJhdGl2YQ0KYGBge3J9DQpmb2xoYSA8LSBpbWFnZV9wbGltYW4oInNldl9sZWFmLmpwZyIsIHBsb3QgPSBUUlVFKQ0KIyBtZWFzdXJlX2Rpc2Vhc2VfaXRlcihmb2xoYSwgdmlld2VyID0gIm1hcHZpZXciKQ0KDQpgYGANCg0KDQojIyBVc2FuZG8gcGFsZXRhcw0KYGBge3J9DQpmb2xoYSA8LSBpbWFnZV9pbXBvcnQoInNveWJlYW5fcnVzdC9zb3lfMS5qcGciLCBwbG90ID0gVFJVRSkNCnBhbHMgPC0gaW1hZ2VfaW1wb3J0KHBhdHRlcm4gPSAic29qYV8iLA0KICAgICAgICAgICAgICAgICAgICAgcGF0aCA9ICJzb3liZWFuX3J1c3QiLA0KICAgICAgICAgICAgICAgICAgICAgcGxvdCA9IFRSVUUsDQogICAgICAgICAgICAgICAgICAgICBuY29sID0gMykNCnNldiA8LSANCiAgbWVhc3VyZV9kaXNlYXNlKGZvbGhhLA0KICAgICAgICAgICAgICAgICAgaW1nX2hlYWx0aHkgPSBwYWxzJHNvamFfaC5wbmcsICAgICAjIHBhbGV0YSBwYXJhIGZvbGhhIHNhZGlhDQogICAgICAgICAgICAgICAgICBpbWdfc3ltcHRvbXMgPSBwYWxzJHNvamFfcy5wbmcsICAgICMgcGFsZXRhIHBhcmEgZm9saGEgZG9lbnRlDQogICAgICAgICAgICAgICAgICBpbWdfYmFja2dyb3VuZCA9IHBhbHMkc29qYV9iLnBuZykgICMgcGFsZXRhIHBhcmEgZnVuZG8NCnNldiRzZXZlcml0eQ0KDQpgYGANCg0KDQojIyBVc2FuZG8gw61uZGljZXMgZGUgaW1hZ2VtDQpgYGB7cn0NCiMgw41uZGljZSBwYXJhIHNlZ21lbnRhciBhIGZvbGhhIGRvIGZ1bmRvLCBkZXBvaXMgYSBkb2Vuw6dhIGRhIGZvbGhhDQppbWFnZV9zZWdtZW50X2l0ZXIoZm9saGEsIA0KICAgICAgICAgICAgICAgICAgIGluZGV4ID0gYygiQiIsICJOR1JESSIpLA0KICAgICAgICAgICAgICAgICAgIG5jb2wgPSAzKQ0KDQojIElNUE9SVEFOVEU6IHRoZXNob2xkIHZhaSBhbHRlcmFyIGEgc2V2ZXJpZGFkZQ0Kc2V2MiA8LSANCiAgbWVhc3VyZV9kaXNlYXNlKGZvbGhhLA0KICAgICAgICAgICAgICAgICAgaW5kZXhfbGIgPSAiQiIsICAgICAgICAgICAgIyDDrW5kaWNlIHBhcmEgc2VnbWVudGFyIGEgZm9saGEgZG8gZnVuZG8NCiAgICAgICAgICAgICAgICAgIGluZGV4X2RoID0gIk5HUkRJIiwgICAgICAgICMgw61uZGljZSBwYXJhIHNlZ21lbnRhciBmb2xoYSBkb2VudGUgZSBzYWRpYQ0KICAgICAgICAgICAgICAgICAgdGhyZXNob2xkID0gYygiT3RzdSIsIDApKSAgIyB0aHJlc2hvbGQgcGFyYSBvcyBkb2lzIMOtbmRpY2VzLCByZXNwZWN0aXZhbWVudGUNCnNldjIkc2V2ZXJpdHkNCg0KIyBlbSBsb3RlICh1dGlsaXphbmRvIGFzIHBhbGV0YXMpDQpzZXYzIDwtIA0KICBtZWFzdXJlX2Rpc2Vhc2UocGF0dGVybiA9ICJzb3lfIiwNCiAgICAgICAgICAgICAgICAgIGRpcl9vcmlnaW5hbCA9ICJzb3liZWFuX3J1c3QiLCAgICAjIHBhc3RhIGNvbSBhcyBpbWFnZW5zIG9yaWdpbmFpcw0KICAgICAgICAgICAgICAgICAgZGlyX3Byb2Nlc3NlZCA9ICJwcm9jX3J1c3QiLCAgICAgICMgcGFzdGEgY29tIGFzIGltYWdlbnMgcHJvY2Vzc2FkYXMNCiAgICAgICAgICAgICAgICAgIHNhdmVfaW1hZ2UgPSBUUlVFLA0KICAgICAgICAgICAgICAgICAgcGxvdCA9IEZBTFNFLA0KICAgICAgICAgICAgICAgICAgaW1nX2hlYWx0aHkgPSBwYWxzJHNvamFfaC5wbmcsDQogICAgICAgICAgICAgICAgICBpbWdfc3ltcHRvbXMgPSBwYWxzJHNvamFfcy5wbmcsDQogICAgICAgICAgICAgICAgICBpbWdfYmFja2dyb3VuZCA9IHBhbHMkc29qYV9iLnBuZywNCiAgICAgICAgICAgICAgICAgIHBhcmFsbGVsID0gVFJVRSkNCiMgc3RhbmRhcmQgYXJlYSBkaWFncmFtcw0Kc2FkKHNldjMsIG4gPSA2KQ0KDQpzZXY0IDwtIA0KICBtZWFzdXJlX2Rpc2Vhc2UocGF0dGVybiA9ICJzb3lfIiwNCiAgICAgICAgICAgICAgICAgIGRpcl9vcmlnaW5hbCA9ICJzb3liZWFuX3J1c3QiLA0KICAgICAgICAgICAgICAgICAgcGxvdCA9IEZBTFNFLA0KICAgICAgICAgICAgICAgICAgaW5kZXhfbGIgPSAiQiIsDQogICAgICAgICAgICAgICAgICBpbmRleF9kaCA9ICJOR1JESSIsDQogICAgICAgICAgICAgICAgICB0aHJlc2hvbGQgPSBjKCJPdHN1IiwgMCksDQogICAgICAgICAgICAgICAgICBwYXJhbGxlbCA9IFRSVUUpDQoNCg0KIyBwYWxldGFzIHZzIMOtbmRpY2VzIGRlIGltYWdlbQ0KDQpkZiA8LSANCiAgZGF0YS5mcmFtZShwYWxldGFzID0gc2V2MyRzZXZlcml0eSRzeW1wdG9tYXRpYywNCiAgICAgICAgICAgICDDrW5kaWNlcyA9IHNldjQkc2V2ZXJpdHkkc3ltcHRvbWF0aWMpDQpkZg0KDQpnZ3Bsb3QoZGYsIGFlcyhwYWxldGFzLCDDrW5kaWNlcykpICsNCiAgZ2VvbV9wb2ludCgpICsNCiAgZ2VvbV9hYmxpbmUoaW50ZXJjZXB0ID0gMCwgc2xvcGUgPSAxLCBjb2xvciA9ICJyZWQiKSArDQogIHRoZW1lX2J3KGJhc2Vfc2l6ZSA9IDE2KSArDQogIHhsaW0oMCwgMTAwKSArDQogIHlsaW0oMCwgMTAwKSArDQogIGNvb3JkX2VxdWFsKCkNCg0KYGBgDQoNCiMjIEVtIGdyaWRzDQoNCmBgYHtyfQ0KZm9saGFfZ3JpZCA8LSBpbWFnZV9pbXBvcnQoImRpc2Vhc2Vfc2hwLmpwZyIsIHBsb3QgPSBUUlVFKQ0KDQpzZXY1IDwtIA0KICBtZWFzdXJlX2Rpc2Vhc2Vfc2hwKGZvbGhhX2dyaWQsDQogICAgICAgICAgICAgICAgICAgICAgaW5kZXhfbGIgPSAiSFVFMiIsDQogICAgICAgICAgICAgICAgICAgICAgaW5kZXhfZGggPSAiTkdSREkiLA0KICAgICAgICAgICAgICAgICAgICAgIHRocmVzaG9sZCA9IGMoIk90c3UiLCAwKSwNCiAgICAgICAgICAgICAgICAgICAgICBmaWx0ZXIgPSAxMCwNCiAgICAgICAgICAgICAgICAgICAgICBuY29sID0gMywNCiAgICAgICAgICAgICAgICAgICAgICBucm93ID0gMywNCiAgICAgICAgICAgICAgICAgICAgICBwbG90ID0gRkFMU0UpDQpwbG90KGZvbGhhX2dyaWQpDQpwbG90KHNldjUkc2hhcGVmaWxlcykNCnBsb3RfbWVhc3VyZXMoc2V2NSwgInN5bXB0b21hdGljIikNCg0KYGBgDQoNCg0KDQo=